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ABSTRACT 

A demonstration project has been undertaken in the 
Pacific Northwest to determine the applicability of 
satellite snowoover observations for operational use 
in three test areas of the Columbia River Basin. 

INTRODUCTION 

Operational application of satellite snowoover observations 
in the Northwestern United States is one of several demonstration 
projects sponsored by the National Aeronautics and Space Admin- 
istration’s Goddard Space Flight Center. A contract was entered 
into between NASA and the Bonneville Power Administration and the 
North Pacific Division, Corps of Engineers. BPA is the primary 
contracting agency with Dr. Mark Meier of the Geological Survey 
acting as the technical field advisor. This contract implements 
a study using satellite data for determining snowoover area which 
is a prime factor in forecasting runoff for flood control and 
power operation. 


APPLICATION 

More than 80 percent of the electric energy produced in 
the Pacific Northwest comes from hydroelectric power generation. 

As of December 31 , 1974, there are 159 hydroelectric plants which 
generate over 21 million kilowatts and 39 thermal plants which 
generate 3.6 million kilowatts. An additional 8.6 million kilo- 
watts of hydro is under construction with another 5.0 million 
authorized or licensed. There are 2.4 million kilowatts of ther- 
mal under construction and 8.6 million kilowatts licensed. Exist- 
ing, under construction, and authorized or licensed hydro plants 
total 34.6 million kilowatts and thermal plants 14.6 million kilo- 
watts. 

The hydro-based system in the Pacific Northwest makes this 
area second to none as a producer of hydroelectric power. The 
hydroelectric projects of the region have been undertaken by pri- 
vate utilities, public agencies, municipalities, and the United 
States Federal Government. These developments take into consid- 
eration all of the multipurpose uses such as power, flood control, 
irrigation, recreation, pollution control, and navigation. 
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Although considerable development has been undertaken by all par- 
ticipating agencies, the larger multipurpose developments have 
been constructed by the Corps of Engineers and the Bureau of Rec- 
lamation. Nearly all of the economically feasible hydro sites 
have been developed and thermal power plants will be needed more 
and more to meet the ever growing power load of the area. Even- 
tually thermal will carry the base load with hydro providing the 
peak power requirements. There are still seme feasible flood con- 
trol sites in the area. Under all of these conditions, an accu- 
rate forecast of voltxne runoff for proper scheduling of the river 
operation to meet all of the multipurpose project uses is of 
paramount importance, Figure 1. 

The Oolunbia River Basin contains more than 668,200 km 2 
(258,000 mi 2 ), almost twice the size of California. The Basin 
contains meet of Washington, Oregon, and Idaho; that part of Mon- 
tana west of the Rocky Mountains; small areas of Wyoming, Utah, 
and Nevada; and the southeastern part of British Oolunbia. 

The Columbia River is second only to the Mississippi in 
average runoff for rivers in the United States. The average 
annual precipitation over the Basin is about 71 cm (28 in) . Of 
this amount, 30.5 cm (12 in) is returned to the atmosphere by 
evapo transpiration and 2.54 cm (1 in) of the amount withdrawn for 
beneficial use is censured leaving 38 cm (15 in) for runoff. This 
amount is equivalent to 228 billion m 3 (185 million acre-feet) of 
water annually or a flow of 7220 m 3 per second (255,000 cubic feet 
per second) . Average flows for key power points in the Basin are 
shown in Table 1. 

Flow records of the Oolunbia River at The Dalles, Oregon, 
cover a period of over 96 years. During this time, the flows have 
ranged from a low of 991.1 m 3 /s (35,000 ft 3 /s) on January 12, 1937, 
to a high of 35,112.9 m 3 /s (1,240,000 ft 3 /s) on June 6, 1894. The 
high flow is about six times the long-term average of 5493.5 ra 3 /s 
(194,000 ft 3 /s) at The Dalles and the low flew is about one-sixth 
the average. 

Storage projects in the Basin sure used to regulate the 
streamflow runoff. There are 25.5 billion m 3 (20.7 million acre- 
feet) of storage in Federal reservoirs, 5.8 billion m 3 (4.7 mil- 
lion acre- feet) in ncn-Federal reservoirs, and 19.1 billion m 3 
(15.5 million acre-feet) in the reservoirs in the Canadian portion 
of the Basin. Thus, usable storage capacity for power projects 
existing, under construction, and authorized or licensed totals 
about 50.4 billion m 3 (41 million acre-feet). Optimal storage in 
the Basin is estimated to be serrewhat in excess of 61.7 billion m 3 
(50 million acre-feet) . The ratio of runoff to storage available 
provides a measure of the usability of the water available. A 
oarparison of the Oolunbia, the Colorado, and the Missouri project 
storage capacities is shown in Table 2. From these figures it can 
be seen that the capacity of the Colorado and Missouri projects 
are several times the average annual runoff. Conversely, the res- 
ervoir storage capacity in the Oolurbia River system is about one- 
quarter the average runoff. The low ratio of reservoir capacity 
to average runoff in the Oolurbia systan can produce large amounts 
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tabu: 1 


Average Flow for Period July 1928, Through June 1968, 
Adjusted to 1970 Level of Development 



Annual Runoff 



Annual Runoff 



Million 



Million 


Project 

Acre-Poet 

kefs 

Project 

Acre-Feet 

kefs 

MOUTH OF 



N. F. d£AFWATER RIVER 


CQLLMJIA RIVER 

185.0 

255.0 


' 





Dworshak 

4.1 

5.6 

Bonneville 

132.8 

183.3 




The Dalles 

128.9 

177.9 

M. F. CIZAJWATER RIVER 


John Day 

124.9 

172.4 




MoNary 

123.0 

169.8 

Penny Cliffs 

4.9 

6.7 

Ben Franklin 

85.9 

118.5 




Priest Rapids 

85.8 

118.4 

GRAND HCNDE RIVER 



Wanapun 

85.7 

118.3 




Rock Island 

85.6 

118.2 

Wenaha 

2.1 

2.9 

Rocky Reach 

83.1 

114.7 




Wells 

81.5 

112.5 

SA1AEN RIVER 



Qiief Joseph 

78.2 

108.0 




Grand Coulee 

78.0 

107.7 

Lower Canyon 

7.8 

10.8 

Murphy Creek 

51.4 

70.9 

Freedom 

7.7 

10.6 

Hugh Keenleyside 

29.0 

40.2 

Crevioe 

7.0 

9.7 

Revelstoke Canyon 

20.9 

28.8 




Downie Creek 

18.7 

25.8 

CHELAN RIVER 



Mica 

14.9 

20.5 




Calamity Curve 

4.7 

6.5 

Chelan 

1.4 

2.0 

Luxor 

1.4 

1.9 

SPOKANE RIVER 



EESCHUIES RIVER 



Long iAke 

5.7 

7.8 

Pelton 

3.1 

4.3 

Nine Mile 

5.2 

7.2 




Mairoe Street 

4.9 

6.7 

SNAKE RIVER 



Post Falls 

4.5 

6.2 

Ice Harbor 

34.6 

47.7 

POD ORIELLE RIVER 



lower Mcnunental 

34.6 

47.7 




Little Goose 

34.2 

47.2 

Waneta 

20.1 

27.8 

Lower Granite 

34.1 

47.1 

BotnxJary 

19.4 

26.8 

Asotin 

23.1 

31.9 

Box Canyon 

18.8 

26.0 

China Gardens 

20.6 

28.5 

Albeni Falls 

18.3 

25.3 

^jpaloosa 

12.4 

17.1 

Cabinet Gorge 

15.8 

21.8 

Oxbow 

12.0 

16.5 

Nooocn Rapids 

14.1 

19.4 




Thcrpeon Falls 

14.3 

19.8 

SNAKE RIVER TRIBUTARIES 





CIZAFWATER RIVER 



FLATHEAD RIVER 



Lenore 

10.2 

14.1 

Buffalo topids No. 

4 8.5 

11.8 




Kerr 

8.0 

11.5 
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TABLE 1 Average Flow for Period July 1928, Through June 1968, 
Adjusted to 1970 Level of Development 
(continued) 


Annual Runoff 



Annual Runoff 



Million 



Million 


Project 

Acre- Feet 

kefs 

Project 

Acre-Feet 

kefs 

S. F. FLATHEAD RIVER 



M3S7JZIE RIVER 



Hungry Horse 

2.5 

3.5 

weiterville 

3.3 

4.5 




Leaburg 

3.1 

4.3 

KocnrmY rivfp 









S. F. M3CEHZIE RIVER 


Brilliant 

22.2 

30.6 




Oorra Linn 

20.0 

27.6 

Cougar 

0.6 

0.8 

Katka 

10.2 

14.1 




Kootenai Falls 

9.0 

12.4 

M. F. WILLAMETTE 

RIVER 


Libby 

8.2 

11.3 




Bull River 

4.4 

6.1 

Lookout Point 

2.1 

2.9 




Hills Creek 

0.8 

1.1 

DUNCAN RIVER 



COWLITZ RIVER 



IXjncan Reservoir 

2.5 

3.5 

Mayfield 

4.4 

6.1 

WILLAMETTE RIVER 



Mossyrock 

3.7 

5.1 

Sullivan 

22.2 

30.6 

LEWIS RIVER 



CLACKAMAS RIVER 



Ariel 

3.5 

4.8 




Yale 

2.8 

3.9 

North Fork 

2.0 

2.7 

Swift No. 2 

2.1 

2.9 

OAK GRCVE FORK 



BAKER RIVER 



Oak Grove 

0.4 

0.5 

lower Baker 

1.9 

2.6 

Timothy Meadows 

0.7 

0.1 

Upper Banker 

1.4 

2.0 

NORTO SANTIAM RIVER 



SKAGIT RIVER 



Detroit 

1.5 

2.1 

Gorge 

3.3 

4.5 




Diablo 

3.0 

4.1 

SOUTO SANTIAM RIVER 



Ross 

2.5 

1.4 

Foster 

1.8 

2.5 

NI SOU ALLY RIVER 



MIDDLE SANTIAM RIVER 



La Grande 

1.0 

1.4 

Green Peter 

1.2 

1.6 

SKOKCMISH RTVER 




Cushman Mol 2 0.6 0.8 


BPA - Branch of Power Resources 
April 22, 1975 


RISPRODUCIBILrrY OF THE 
ORIGINAL PAGE It ?O0* 
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TABI£ 2 Flew and Storage - Columbia, Colorado and Missouri Rivers 



Storage 
Capacity 
Billion m 3 1/ 


Annual 

Runoff 

Billion m 3 1/ 



Headwater 
Elevation 
meters 2/ 

Drainage 

Area 

Billion m 2 3/ 



Average 

Maximum 

Minimum 



Columbia 

50.6 (41) 

228.2 (185) 

315.8 (256) 

149.3 

(121) 

807.7 (2.65) 

668.2 (258) 

Colorado 

66.6 (54) 

17.3 (14) 

27.1 (22) 

4.9 

(4) 

2652.0 (8.70) 

626.8 (242) 

Missouri 

92.5 (75) 

29.6 (24) 

45.6 (37) 

13.6 

(11) 

2134.0 (7.00) 

1364.9 (527) 


1/ Figures in ( ) are x 10 6 acre-feet. 

7/ Figures in ( ) are mean sea level in feet x 10 3 . 

1 / Figures in ( ) cure miles 2 x 10 3 . 



of secondary energy. This low storage to runoff ratio and the 
large difference in average and maximum streamflow at The Dalles, 
Oregon, shows the need for accurate volume runoff forecasting and 
optimum use of storage for flood regulation. 

The National Weather Service (NWS) , the Corps of Engineers, 
North Pacific Division (USCE, NPD) , and the Bonneville Power 
Administration (BPA) are cooperatives in the Columbia River Fore- 
casting Service (CRFS) . The overall goal of the CRFS is to pool 
certain resources of the agencies in the interest of improving 
streamflow forecasting methods, to provide uniform forecasts, and 
to increase the efficiency of operation. 

Daily operational runoff forecasts for streams in the 
Pacific Northwest are made unsing the Streamflow Synthesis and 
Reservoir Regulation (SSARR) computer model. The program descrip- 
tion and user manual along with associated publications referenced 
in the manual provide detailed background on the development of 
the model and hew it is used for operational river forecasting and 
river management activities (1) . 

Cne of the major factors in runoff in the Columbia River 
Basin is snowmelt. Snowmelt calculation with the SSARR model is 
made by the temperature index method or by the use of the gener- 
alized snownelt equation for a partly forested area. At the pres- 
ent time we are using the temperature index method for day-to-day 
forecasts. In general, the snowmelt equation is not used for 
daily operational forecasts because of the lack of real time 
energy budget data. 

Snowmelt runoff in inches is determined for the snow area 
using the temperature index method: 

”■ < t a-V r [§r] 

where 

m = Snowmelt runoff in inches of water over the -snowoover 
area. 

T a * Period tsrperature at the median elevation of the melt- 
ing snowpack. 

= Base temperature (°F) , specified as a constant for a 
watershed. 

R = Melt rate, specified to the oemputer, or given as a 
function of accumlated runoff, in inches of water per 
degree day. 

PH = Period length in hours. 

Values used for the base temperature and the melt rates can be 
adjusted to either minimum, mean, or maximum daily temperature 
daily temperature data. Melt rates can be specified for each day 
of the run to conform to the natural variability encountered dur- 
ing the melt season. 

TVoo methods are available to evaluate snowmelt from a 
watershed. The basin method evaluates the snow-covered area run- 
off relationships using a snowoover depletion function. The other 
method divides the watershed into elevation bands with each band 
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being examined separately with respect to snow accumulation and 
melt, studies are now underway to evaluate this second method 
for some of the sub-basins in the region. 

Snowoover observations in the United States portion of the 
Basin are made by the Corps of Engineers personnel in snail air- 
craft flying at a low altitude. Similar flights are made in the 
Canadian portion of the Basin by the British Columbia Hydro and 
Power Authority. Observations of this type are made by experi- 
enced personnel and must, of necessity, be subjective. Two to 
four flights are generally made each season for each area depend- 
ing upon the type of season and flying conditions. In the United 
States portion of the Basin, flights generally occur in April, 

May, and June and in the Canadian portion in May and June. Though 
quite effective, these flights may present a safety problem and 
are relatively expensive. 

The satellite data collection progr a m provided an excellent 
opportunity for testing aerial snowcover observations as an alter- 
nate to such determination by small aircraft. BPA and the Corps 
of Engineers formulated a proposal to the National Aeronautics and 
Space Mministration/Goddard Space Flight Center (NASA/GSPC) to 
study snowoover observations which would hopefully improve fore- 
casting procedures. 

The primary contracting agency was BPA. A Memorandum of 
Understanding was signed to cover contractual work for one year by 
BPA and the Corps of Engineers for NASA/GSFC. Work included: (1) 

the mapping of snowoover in three candidate areas for the 1973 and 
1974 spring melt seasons, (2) input of snowoover information in 
the SSARR model for reconstitution of the temporal distribution of 
the spring melt runoff, and (3) compare results with reconstitu- 
tions made with the conventional data for the past two seasons. 

The basins studied are: the Boise River above Lucky Peak 

Dam, drainage area 7254 km 2 (2800 mi 2 ) ; North Santiam River above 
Detroit Dam, drainage area 1132 Jan 2 (437 mi 2 ); and the upper Snake 
River above Palisades Dam, drainage area 13,489 km 2 (5208 mi 2 ). 
These three basins were chosen initially because of some associ- 
ated work being done by Dr. Mark Meier of the Geological Survey on 
the Santiam Basin, reconstitution work already done by the Corps 
of Engineers and the National Weather Service on the Boise and 
Snake River Basins, and because it was expected that adequate sat- 
ellite imagery would be available. Satellite data for 1975 is 
being examined new to determine adequate cloud-free coverage over 
the Flathead River above Hungry Horse Dam, the Kootenai River 
above Libby Dam, and the North Pork of the Clearwater River above 
Dworshak Dam, Figure 2. 

Imagery from the Earth Resources Satellite, LANDSAT-1 , was 
assembled for the 1972-73 and 1973-74 snowmelt season. This 
imagery was processed under a subcontract with Stanford Research 
Institute (SRI) utilizing its Electronic Image Analysis Console 
(ESIAC) to determine the percentage of snowoover area (2). The 
SRI ESIAC greatly expands the capability of the human analysts, 
particularly their ability to effectively use time-sequence, 
multispectral imagery and provide objective and repeatable mea- 
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FIGURE 2 

LANDSAT Orbital Tracks Over Tsst Basins 
Using Canadian Numbering System 



sursnents. It also provides a bridge between manual photo inter- 
pretation and machine processing. SRI presented three main 
options, in ascending order of cost and degree of refinement, for 
analyzing the imagery: manual photo interpretation, ESIAC clas- 

sification from one spectral band, and advanced ESIAC analysis 
using two spectral bands and use of variable area, above thresh- 
old elevation masks to establish localized snowline elevations and 
to predict shadow regions. The second option using ESIAC clas- 
sification from one spectral band (MSS- 5) was chosen. 

The satellite orbited tracks using the Canadian numbering 
system and the frame coverage is shown on Figure 2. One LANDSAT 
frame will generally cover the North Santiam Basin, and this basin 
was used to develop detailed procedures. About two frames are 
needed to cover the Boise Basin and from two to four frames are 
required for the Upper Snake Basin. 

Snowoover area was determined from a program that combines 
an objective classification using single-band (MSS-5) radiance 
thresholding with subjective editing using elevation contours and 
other reference data as a guide. LANDSAT imagery was scaled to 
the ESIAC TV screen, about 55 x 40 km or 0.5° by 0.5° longitude 
and latitude. This scaling approximately matched the display 
resolution to the LANDSAT imagery resolution. This magnification 
required one TV view for the North Santiam Basin, four views for 
the Boise Basin and eleven views for the Upper Snake Basin. A 
binary basin outline nap and a Universal Transverse Mercator (IHM) 
grid in 2.5 km x 2.5 km boxes were superimposed over each view. 

All of these data were spacially registered to each other during 
the entry process. 

Registered color, time-lapse sequences were created for 
each basin. A mask of the snowoover area was generated for each 
view based on radiance thresholding, localized threshold adjust- 
ments, and superimposed elevation contours with final editing when 
necessary. The masks were documented by photography, by numeri- 
cal pixel count of the total area, and by a digital array depict- 
ing tenths of snowoover for each 2.5 x 2.5 km cell. The cell-by- 
oell numerical documentation provides a convenient means for 
machine recording and manipulating the data, comparing with other 
methods of determining snowoover, and, more importantly, for use 
with advanced modeling techniques such as elevation bands as men- 
tioned earlier. 

Fined, results of the work on the measurement of snowoover 
from satellite imagery are not as yet available. Very adequate 
data were obtained for the North Santiam Basin; however, this 
basin was chosen primarily for the work being done by Dr. Marie 
Meier to test snowoover area methods. The North Santiam Basin is 
a coastal basin and the rmoff is, for the most part, during the 
winter season from rain. Thus, it is not included as a primary 
snow runoff basin in the SSARR model. Data for the North Santiam 
Basin is shown in Table 3. 

The Upper Snake Basin snowoover is shown in Table 4. Dur- 
ing 1973, because of persistent cloud cover, LANDSAT- 1 images were 
usable for only three dates, 29 March, 22 May, and 9 June. More- 
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TABLE 3 N. Santiam Basin Sncw Coverage 


Basin Area = 1146 km 2 Average Elevation 3900 ft MSL 


Date 

% of Basin Area 
Covered by 

Snow Coverage 

Usable imagery 

Area (km 2 ) 

% of Basin Area* 

1973 




6 January 

100 

858.4 

74.9 

11 February 

100 

891.4 

77.8 

6 April 

100 

519.1 

45.3 

24 April 

100 

475.6 

41.5 

12 May 

100 

205.1 

17.9 

1974 




1 January 

100 

975.2 

85.1 

6 February 

100 

1002.8 

87.5 

24 February 

100 

1080.7 

94.3 

12 June 

100 

405.7 

35.4 

30 June 

100 

190.2 

16.6 


♦Derived by dividing the measured snow area within the 
basin by the viewable area within the basin. 
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TABLE 4 Upper Snake Basin Snow Coverage 
Basin Area = 13,339.6 km 2 Average Elevation ^ 8000 ft MSL 


Date 


Snc* 

Area (km 2 ) 

r Coverage 
% of Basin Area* 




LANDSAT 

SSAKR 

1973 





29 March 

96.9 

11888.0 

92.0 

100 

22 May 

95.0 

5128.4 

40.5 

33 

9 June 

94.9 

1738.8 

13.7 

19 

1974 





22 June 

96.8 

1942.2 

15.3 



♦Derived by dividing the measured snow area within the basin 
by the viewable area within the basin. 


82 








over, no aircraft snow flights were made in 1973 as this was a 
very low year for precipitation and runoff. Agreement between the 
SSARR depleted data using the temperature index and the LANDSAT 
data are fair. The large gap between 29 March and 22 May misses 
the high runoff which started in mid-May and peaked on 22 May. In 
1974 only one usable image was available for 22 June. This was 
much too late in the season and efforts are being made to secure 
additional coverage from overlapping imagery to the east that 
would cover most of the basin. 

The Boise Basin snowoover is shown in Table 5. Again, 
there were no aircraft snow flights in 1973 but the SSARR depleted 
coverage shows good agreement with LANDSAT-1 data particularly 
during the runoff period. Of particular interest are the data for 

1974. Very good agreement is evident between the satellite and 
SSARR data except for 14 April where the SSARR data shows 57 per- 
cent and the satellite data 74.1 percent coverage. The probable 
cause of this disagreement was a thin skin of new snow that fell 
cn 12 April. This new snow extended below the 4000 foot elevation 
and showed up cn the satellite; but, because of its transitory 
nature, provided little in the way of runoff. 

Comparison of the actual hydrograph and the reconstitution 
of satellite and SSARR data for runoff from the Boise Basin is 
shown on Figure 3. As can be seen, the correlation is very good. 

It is still too early to draw any definite conclusions on 
the use of satellite data for determing area of snowoover in this 
region; particularly with the large number of cloud-covered days, 
the problem of determining the snowline in forested areas, and the 
lack of high resolution with the Synchronous Meteorological Sat- 
ellite (SMS) data. The addition of new satellites this past year 
and future improvement in resolution should materially improve 
satellite data. Yet to be tested is the turn around time in get- 
ting the satellite data from the receiver sites to the user for 
real-time application. Nonetheless, what we have seen ‘to date 
shows that satellite determination of snow coverage is a promising 
adjunct to our current methods. Continuation of our project, par- 
ticularly with the deep snowpack and late runoff experienced in 

1975, should provide us with much added information. We are much 
indebted to NASA/GSFC for supporting this valuable project. 
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TABLE 5 Boise Basin Snow Coverage 


Basin Area = 7254 km 2 Average Elevation v 6000 ft MSL 


Date 

T mb’' , "' jW 

jaai»,T 

Snow Coverage 

Area (km 2 ) 

% of Basj 

m Area* 




LANDSAT 

SSARR 

1973 





19 April 

100 

4750.4 

65.5 

54 

7 May 

100 

2270.4 

31.3 

32 

12 June 

100 

594.0 

8.2 

6 

30 June 

100 

278.6 

3.7 

2 

1974 





14 April 

100 

5377.0 

74.1 

57 

2 May 

100 

2609.5 

36.0 

34 

25 June 

100 

504.5 

7.0 

6 


♦Derived by dividing the measured snow area within the basin 
by the viewable area within the basin. 
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Figure 3 : Boise Basin above Lucky Peak Dam, Hydrograph and 
SSARR Computed Flows 
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